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МОДЕЛИРОВАНИЕ ПРОЦЕССА ОХЛАЖДЕНИЯ ЭЛЕМЕНТОВ 
РАДИОЭЛЕКТРОННОИ АППАРАТУРЫ В 5ОМОМ/ОВК$ 


С увеличением интеграции современной радиоэлектронной аппаратуры (РЭА) существенно возрастает и тру- 
доемкость ее проектирования, что связано с ростом размерности задач, решаемых в процессе разработки. 
Эффективным способом решения этой проблемы является применение системы автоматизированного проек- 
тирования 5о/АИ/огк$. В работе с помощью этой системы решается задача расчета параметров радиатора 
охлаждения для бесперебойной работы одного из элементов РЭА. 

Ключевые слова: моделирование процесса охлаждения, тепловыделяющие элементы, теплоотдача, метод 
конечных элементов (МКЭ), $оаИ/огк$ Нои/ 5ти/айоп. 


Введение. Блок обработки приемного сигнала (далее блок) с радиатором жидкостного охлажде- 
ния его тепловыделяющих элементов входит в состав бортовой аппаратуры летательных аппара- 
тов. К этой группе продукции предъявляются особые требования по надежности. Надежность, 
прежде всего, определяется качеством разработки конструкции и использованием максимальных 
возможностей средств по моделированию и конструированию. Задачей конструктора блока явля- 
ется оптимизация конструкции с целью обеспечение неизменности параметров измерения при 
дестабилизирующих климатических воздействиях, интенсивность которых связана с условиями 
эксплуатации, при минимизации материалоемкости и сокращении номенклатуры элементов и ма- 
териалов. В частности, температура кристалла интегральной микросхемы не должна превышать 
+70 °С; при этом атмосферное давление - 150—450 мм рт. ст., температура окружающего воз- 
духа —50 9С — +65 9С, влажность воздуха - 98 %. 

Для РЭА разработана методика теплового расчета коэффициентным способом [1, 2], кото- 

рая с достаточной точностью просто и быстро позволяет решать конструкторские вопросы тепло- 
обеспечения радиоаппаратуры. 5оП!а\М\/огК$ — одна из систем, инструменты для моделирования ме- 
ханической части РЭА которой дополняются специальными приложениями, позволяющими прово- 
дить частотный анализ, тепловые и прочностные расчеты. У тепловых расчетов в Зо!а\М/огК$ есть 
свои преимущества: наглядность результатов, автоматическое создание отчета анализа, удобство 
в использовании, исключение ручных расчетов, облегчение интерпретации результатов анализов. 
Цель исследования: моделирование процесса охлаждения блока в 5о!а\М\югк$, оценка эффек- 
тивности охлаждения тепловыделяющих элементов радиатором жидкостного охлаждения в опре- 
деленных климатических условиях; логико-математический выбор устойчивых компромиссов для 
удовлетворения противоречивых требований технического задания по назначению и надежности 
исследуемого объекта. 
Постановка задачи. В состав блока входят тепловыделяющие элементы, нагреваемые за счет 
преобразования электрической энергии в тепловую и повышенной температуры окружающего 
воздуха, корпус блока, радиатор, отводящий тепло от тепловыделяющих элементов на корпус 
блока. К корпусу блока примыкает радиатор с циркулирующей охлаждающей жидкостью. 

В расчетном модуле бойа\Мотщ$ Ном/ УтиаНоп необходимо задать условия задачи, макси- 
мально приближенные к реальным условиям эксплуатации, провести расчеты и интерпретацию 
полученных результатов. 

В основе расчетной методики лежит, во-первых, математическая модель рассчитываемых 
физических процессов и, во-вторых, способ решения поставленной математической задачи. 

Математические модели конкретной инженерной задачи представляют собой начально- 
краевые задачи, отражающих законы физики и механики с привлечением полуэмпирических 
и эмпирических констант и зависимостей. 

Движение и теплообмен в жидкости и газе моделируется с помощью уравнений Навье- 


Стокса [31, описывающих законы сохранения массы, импульса и энергии этих сред. Эта система 
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уравнений нестационарного пространственного течения имеет следующий вид в рамках подхода 
Эйлера в декартовой системе координат (х, г = /, 2, 3): 
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где { — время; и — скорость среды; р— плотность среды; Р’— давление среды; 5,— внешние 
массовые силы (действие гравитации), действующие на единичную массу среды; Ё— пол- 
ная энергия единичной массы среды; О, — тепло, выделяемое тепловым источником 


в единичном объеме среды; т, — тензор вязких сдвиговых напряжений; 4,— диффузион- 


И 
ный тепловой поток. 
Уравнение состояния газообразной среды: 


р=р(Р,Т, у} (4) 
где у=(у,у,, ..., у, } — вектор концентраций компонентов среды. 
Уравнение СОСТОЯНИЯ ЖИДКОСТИ: 
р=Р(Т, у). (5) 
Тензор вязких сдвиговых напряжений определяется следующим образом: 
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где и=н, +н., — коэффициент динамической вязкости и; и,— коэффициент турбулентной 
вязкости; д, — дельта-функция Кронекера; А — кинетическая энергия турбулентности. 
Диффузионный тепловой поток моделируется с помощью уравнения 
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где с, =0,9, Р.— число Прандтля; с›— удельная теплоемкость при постоянном давлении; 


а 
Т - температура среды. 
Теплопередача в твердых телах рассчитывается с помощью уравнения [4]: 


дрет _ д ег 0. г 





[6] дх. | дх. 


р 
где с — удельная теплоемкость; 1’ — температура; ^. — теплопроводность. 
Радиационный теплообмен моделируется в соответствии с законом Стефана—Больцмана: 
Он = оо (14 и (9) 
где = — степень черноты поверхности; с, — постоянная Стефана—Больцмана; 7’, — температура 


поверхности, испускающей излучение; Г, — температура окружающей среды. 
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Поскольку используемые в математической модели системы дифференциальных уравне- 

ний не имеют аналитического решения, они приводятся к дискретному виду и решаются на неко- 
торой расчетной сетке. 
Моделирование в $ойа\МогКк$. В системе 5ойА\М\огк$ (Договор № И05110-77 от 05.11.2009 
г.дстроим объемную модель корпуса и радиатора (рис.1). Для проведения расчета модель блока 
упрощаем в силу лимитированности вычислительных ресурсов: заменяем элементы сложной гео- 
метрии более простыми, исключаем элементы, влиянием на результаты расчета которых можно 
пренебречь. Выполняем две заглушки в радиаторе жидкостного радиатора, «герметизирующие» 
полость каналов внутри радиатора, и две заглушки на корпусе блока. Это сделано для того, что- 
бы сформировать подобласть, занятую охлаждающей жидкостью, в первом случае, и подобласть, 
заполненную воздухом, во втором случае. Условия входа и выхода теплоносителей должны быть 
приложены к внутренним сторонам заглушек. 





а) 6) 


Рис.1. Блок обработки приемного сигнала: 
а- вид спереди; 6 - вид сзади (на виде спереди крышка скрыта) 


Конструкция блока (рис.2) включает в себя корпус 3, тепловыделяющие элементы, креп- 
ления, пластинчатые радиаторы между корпусом и тепловыделяющими элементами. Тепловыде- 
ляющими элементами являются микросхемы. Две микросхемы 1 входят в состав модуля электро- 
питания \УТСОВ. \24С5Т50ВЕЗ с тепловыделением 7,5 Вт каждая и располагаются между платой 
модуля и пластинчатым радиатором. Четыре микросхемы 2, мощностью 25 Вт каждая, входят 
в состав устройства управления и цифровой обработки. Вся конструкция изготовлена преимуще- 
ственно из алюминиевого сплава АМг5, за исключением плат из стеклотекстолита. Радиатор ох- 
лаждения примыкает ко дну корпуса блока. Блок проектируем таким образом, чтобы тепловыде- 
ляющие элементы плотно прилегали к пластинчатым радиаторам, отводящим тепло к корпусу 
блока. А затем от корпуса блока тепло передается общему радиатору жидкостного охлаждения. 
Жидкость, в качестве которой выбран водный раствор этиленгликоля ОЖ-65, циркулирует по 
П-образному каналу круглого сечения диаметром 10 мм. Температура охлаждающей жидкости на 
входе в канал равна +35-5 9С, расход - 1 л/мин. Задачей охлаждения является обеспечение от- 
вода тепла с теплонагруженных элементов на несущую конструкцию, а затем осуществление кон- 
вективной теплопередачи между конструкцией и теплоносителем. 
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Рис.2. Модель блока для расчета: 
а- вид спереди; 6 - вид сзади (на виде спереди крышка скрыта) 


Деградационные процессы при пониженной температуре, связанные с конденсацией вла- 
ги, изменением электрических характеристик, деформаций сопряженных деталей, в расчетах не 
рассматриваются вследствие ограниченных возможностей программы. В расчете рассмотрен слу- 
чай, соответствующий наибольшей вероятности появления перегрева, при котором температура 
окружающего воздуха +65 9С. 

Для адаптации математической модели к конкретной физической задаче и области про- 
странства задаем условия задачи. 

Граничные условия: температура блока — +65 9С, температура окружающего воздуха 
— +65 9С. Назначены две текучие среды: внутри блока: воздух и теплоноситель ОЖ-65 (в каналах 
радиатора). Условие входа охлаждающей жидкости в канал радиатора: температура жидкости 
— +40 °С и скорость жидкости - 0,213 м/с. Условие выхода охлаждающей жидкости: статическое 
давление. На стенках блока — теплоотдача наружных стенок блока с учетом окружающей темпе- 
ратуры, статическое давление - на боковой стенке блока (имитация негерметичности блока). 

Назначены коэффициенты излучения и степень черноты элементов конструкции, массо- 
вый расход, шероховатость поверхности каналов радиатора, мощность термоактивных элементов 
по условию задачи. Создана модель контактного теплового сопротивления между термоактивны- 
ми элементами и корпусом блока, соответствующего защитному покрытию корпуса блока, приме- 
няемому в связи с повышенной влажностью окружающей среды. Также активизирована опция 
«Давление в точке отсчета», при которой статическое давление жидкой среды автоматически 
корректируется по высоте с учетом гравитационного ускорения. 

Для нахождения искомого численного решения задачи континуальная математическая мо- 
дель физических процессов дискретизируется по пространству. Решение задачи о движении 
и теплообмене жидкой среды и газа определяется как установившееся во времени. Исходя из ус- 
ловия Куранта [5], шаги дискретизации по пространству в ячейке определяются в зависимости от 
значений физических величин. На сформированной таким образом некоторой расчетной сетке 
дискретное решение поставленной дифференциальной математической задачи в общем случае 
зависит от размеров ячеек расчетной сетки. Для дискретизации дифференциальных уравнений 
в 5ойА\М/огк$ Ном Ути!аНоп используется метод конечных объемов. Чтобы решить поставленную 





469 


Физико-математические науки 








математическую задачу с заданной точностью, необходимо провести несколько расчетов на раз- 
ных более редких и более частых расчетных сетках, чтобы определить такую частоту сетки, на- 
чиная с которой решение задачи перестает значимо зависеть от частоты сетки, что указывает на 
достижение необходимой точности решения математической задачи, т.е. сеточной сходимости 
решения математической задачи. Таким образом, допустимый уровень разрешения сеткой высту- 
пов поверхности задаем максимальным. Включаем опцию дробления сетки в узких каналах и ус- 
танавливаем критическое число ячеек поперек узкого канала, равное 11. Минимальный размер 
ячеек - 0,123 мм, минимальная толщина стенок твердого тела- 0,7 мм. 
Численные результаты. Решение системы уравнений (1)-(3) для стационарных процессов за- 
висит от следующих условий: 

— физические условия, характеризующие теплофизические свойства среды: 


Сыр», Ву н,; 

_ геометрические условия, определяющие форму и размеры системы, в которой протекает 
процесс; 

— граничные условия, определяющие взаимодействие между различными средами на 
внешних и внутренних границах рассматриваемой системы. 

Варьируя эти условия, можно повлиять на параметры теплового процесса. 

Расчет показал, что при пониженном атмосферном давлении и повышенной температуре 
окружающего воздуха излучение тепла с теплонагруженных элементов в воздух по сравнению 
с нормальными условиями эксплуатации снижается на 17,8 %, а конвекция с воздухом - на 9 %, 
так как уменьшается теплоемкость воздуха. Таким образом, наибольший отвод тепла осуществлен 
радиатором жидкостного охлаждения. Температура микросхем не превышает допустимую и равна 
+52 °С (рис.3), что ниже допустимого значения на 18 9С; температура плат - 54 °С; температура 
воздуха внутри блока - 63,8 °С (рис.4,а); температура жидкости на выходе из радиатора повыси- 
лась на 0,4 9С (рис.4,6). Это позволяет прогнозировать надежную работу блока. Эффективное 
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Рис.3. Распределение температуры в блоке: 
а- вид спереди; 6 - вид сзади (крышка и платы скрыты) 
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охлаждение обусловлено оптимальным сочетанием следующих свойств охлаждающей жидкости, 
величины которых входят в уравнения (1)-(6): коэффициент динамической вязкости, коэффици- 
ент турбулентной вязкости, скоростной напор, плотность, теплоемкость. 
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а) 
Рис.4. Температура воздуха в блоке и потока жидкости в радиаторе охлаждения: 
а- вид спереди; 6 - вид сзади (радиатор скрыт для визуализации результатов) 
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Если варьирование физическими и граничными условиями невозможно, можно изменить, 
во-первых, геометрию канала и, во-вторых, поверхность съема тепла. 

В первом случае при сужении диаметра канала от 11 до 2 мм происходит интенсификация 
теплоотдачи, обусловленная повышением турбулизации пограничного слоя охлаждающей жидко- 
сти вследствие увеличения скорости ее течения. При диаметре меньше 2 мм происходит умень- 
шение теплоотдачи, вызванное сокращением площади съема тепла (см. табл.). При увеличении 
диаметра канала более 13 мм теплоотдача снижается в результате снижения скорости течения. 
При любом из диаметров канала рассматриваемого диапазона теплоотдача эффективна. 


Зависимость параметров теплового процесса от диаметра канала 
































Диаметр о Минимальная Максимальная Скорость течения Отведенная 
канала на выходе температура температура охлаждающей радиатором 
радиатора, мм из канала, °С кристалла, °С кристалла, °С жидкости, М/С тепловая мощность, Вт 

2 40,2 48,5 55,4 5,3 82,1 

4 40,1 48,0 54,5 1,327 81,29 

6 40,15 48,6 55,1 0,589 79,92 

8 40,25 51,3 55,9 0,332 78,02. 

10 40,35 51,7 56,9 0,212 76,39 

11 40,4 52,1 57 0,175 75,38 

13 40,3 53,7 57,3 0,126 76,57 

15 40,45 49,4 55,9 0,094 76,51 




















Во втором случае на теплоотдачу влияет шероховатость поверхности каналов радиатора. 
Причина в том, что вблизи шероховатой стенки скорость потока уменьшается больше, чем вблизи 
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гладкой, что приводит к уменьшению теплового потока на 10 %. Наличие же шероховатости, по 
крайней мере, в модели 5о!а\М\/огк$ Ном ЭитШайоп, никак не сказывается на излучающих свойст- 
вах поверхностей. 
Выводы. Результат расчета показывает наличие в блоке дополнительных напряжений термиче- 
ского происхождения. Величина максимального напряжения составила 1972 Па. Термостойкость 
должны повышать все факторы, уменьшающие величину деформации при термических нагрузках, 
в первую очередь, снижение коэффициента термического расширения и увеличение теплопро- 
водности. Следовательно, для снижения напряжений в конструкции необходимо увеличить тол- 
щину деталей, а также заменить материал конструкции на более термостойкий, например, АМгб. 
Перспектива продолжения работы состоит в решении вопроса о взаимном действии меха- 
нических и тепловых эксплуатационных нагрузок и подтверждении результатов расчетов экспе- 
риментальной проверкой. 
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ЗТМОГАТТОМ ОЕ ЕТЕСТВОМТС ЕООТРМЕМТ ЕТЕМЕМТ$ СООММС 1М $ООМ/ОВК$ 


Напа т Папа мир {Ме едгаНоп тсгеазе оЁ {Пе ир-бю-дае га@ю еесготс едиртепЕ (ВЕЕ) Пеге дго\мз 
эзир$апба!у К5 аедп сотрехКу, СИ 15 аззодае м ап тсгеазе т {Пе айтепзюп оЁ 1е ргоШетз зо№еа 
т {Ле деуеюртепЕ ргосез$. ТНе аррйсаНоп оЁ ЗАВ Зойа\М/огК$ 15 ап еНесёуе мау {© 5о№е {5 ргоМет. Тпе 
ргоЫет оЁ сасща{!та Не рагатеег$ ог Ше гаФаког соо!тд Гог сопбпиоц$ орегаНоп оР опе оГ {Пе еетеп 
ОЕ ВЕЕ 15 эомМеа мВ {Пе пер ог {АБ зу%ет. 

Кеу мога$: соо!та эитшайоп, Гие! еетег(5, Неа! {гапзГег, ИпКе еетепЕ теоа, бо!а\М\/огК$ Ном/ Эит\ЦаНоп. 
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